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INTRODUCTIOX  AND  SUMMARY 


Activity  during  tliis  quarter  included  both  program  develop- 
ment and  additional  data  processing.  The  Andreanof  Islands 
event  (ll/22/()5)  was  processed  througli  the  seismic  source  depth 
determination  program,  with  a scries  of  runs  being  made  using 
various  combinations  of  tlic  14  stations  available.  Although 
four  different  dejUli  ])lot  peaks  were  observed  consistently 
during  these  runs,  further  analysis  eliminated  three  of  them, 
resulting  in  a dei)th  estimate  of  25  km  for  the  event. 

In  the  area  of  program  development,  tw'O  techniques  for 
eliminating  cepstrum  peaks  due  to  the  structure  of  the  source 
spectrum  were  tried.  Both  methods  attempted  to  flatten  the 
source  spectrum  by  inverse  filtering,  but  neither  was  success- 
ful. Tinally,  coding  was  started  on  a new  version  of  the  depth 
determination  program  to  be  implemented  at  tlie  Seismic  Data 
Analysis  Center  (SDAC)  . In  addition  to  allOAving  easier  access 
to  larger  amounts  of  data,  the  SDAC  system  will  provide  greatly 
increased  output  flc.x  ibi  1 ity , designed  to  simjilify  the  analyst's 
task  of  deciding  on  tlie  correct  source  depth. 
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Anch’ca  nof  I si  and  Fvcn  t 

Data  for  the  Andreanof  Island  event  (11/22/65)  consists 
of  seismograms  from  14  stations,  with  A's  ranging  from  34°  to 
73°;  these  seismograms  are  shown  in  Figures  1 and  2.  Recaiise 
of  the  large  number  of  stations  available,  the  depth  determin- 
ation program  could  be  run  using  several  different  groups  of 
stations,  and  the  resulting  depth  plots  could  he  compared  to 
aid  in  deciding  on  the  correct  event  depth.  This  tcchniciue 
prov-^ed  to  he  quite  useful,  since  several  of  the  data  groups, 
including  tlie  one  having  tlie  largest  number  of  stations  (10), 
yielded  four  distinct  depth  plot  peaks.  The  avaailabi  1 ity  of 
results  from  additional  station  groupings  made  it  possilile  to 
resolve  this  ambiguity  and  come  up  with  a single  deptii  esti- 
mate for  the  event. 

Twelve  runs  of  the  deptli  determination  program  were  made 
using  different  station  groupings  and  input  parameters.  Com- 
posite depth  plots  from  tlie  most  informative  of  these  runs 
are  shown  in  Figures  3-18;  results  for  cepstrum  computation 
windows  of  25.6  and  51.2  seconds  arc  shown  for  each  data  set. 
Characteristics  of  the  station  groupings  used  to  obtain  these 
depth  plots  are: 

1.  5 stations,  distributed  over  entire  available 
range  of  A’s. 

2.  10  stations,  distributed  over  entire  available 
range  of  A's. 
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Figure  1 

Andreanof  Islands  Event  Seismograms 
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Composite  Depth  Plot,  Data  Set 


CEPSTRUM  WINDOW  LENGTH  = 51.2  SEC 
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Composite  Depth  Plot,  Data  Set 
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Composite  Depth  Plot,  Data  Set 
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Figure  7 

Composite  Depth  Plot,  Data  Set  2,  Log  Whitening 
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Composite  Depth  Plot,  Data 
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Composite  Depth  Plot,  Data  Set 
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Composite  Depth  Plot,  Data  Set 
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Composite  Depth  Plot,  Data  Set 
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Composite  Depth  Plot,  Data  Set 


16 


r 


TT 


t N 

E 


X 

H 

D, 

Cjp  C'p  rn 

LC  o Q 

CTi  oo 


U 

UJ 

C/2 


VO 

u 

• 

w 

25 

c/: 

<ct 

u 

:r: 

o 

H 

1^ 

w 

II 

C-. 

LD 

o 

>Z 

a 

w 

►— < 

.-J 

<- 

< 

C 

—} 

00 

< 

Cu 

OJ 

O 

O 

H 

i 


LO 

V.  l-UE-O^'.  ClCQO 

o xoa<c/20c/2x;p% 

»— ( I I I t I I I I 

H a,  cj  CL,  cj  ‘.<i  V'.  csi  x; 
< x:CL,L-30CDi;^c_)m 
H 
C/2 


17 


d 


Composite  Depth  Plot,  Data  Set 
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CEPSTRUM  WINDOW  LENGTH  = 25.6  SEC 


7T 


U) 

2:  H < > H oa 

O 2 S 2.  3 2; 

»-H  I I I I I 

H a,  > 2:  2;  0; 

< 2;  E 2:  u;  u 

H 

c/3 


19 


L 


Composite  Depth  Plot,  Data  Set 


B 


o\o 

LO 

Ov 


o\o 

O 

00 


X 

H 

Cu 


Q 


U 

W 

C/D 

CM  U 
• UJ 

rH  00 

to 

II 

C^J 

X o 
{— ' 

O 

II 

w ^ 
t- 

& c: 

o 

C K 
?;  ^ 
f— ( 

s < 

H 

PS  •< 
PD  a 
os 

f-  j 
c/i  < 
CX  H 
W O 
u H 


H < > H CO 

o S'.  PS  X 3 >; 

HH  I I I I I 

E-  D.  > X.  X o; 

< X pp  PS  uj  u 

H 
lO 


20 


Composite  Depth  Plot,  Data  Set 


3. 


6 stations,  A's  within  5°  of  A = 48'’  (middle 
range  A ' s)  . 


4. 

4 stations, 

A's  Av'ithin  5°  of  A 

= 61°  (long 

range  A ' s) . 

5. 

4 stations, 

A's  within  5°  of  A 

= 39°  (close 

range  A's). 

6. 

8 stations. 

arrival  visible  7 

seconds  after  P. 

7. 

5 stations, 

arrival  visible  10 

seconds  after  P 

All  of  the  runs  illustrated  were  done  using  the  first  102.4 
seconds  of  each  seismogram. 

It  is  evident  from  a quick  examination  of  the  plots  in 
Figures  3-18  that  tlie  correct  source  depth  for  this  event  is 
far  from  obvious.  I'our  different  peaks  reajipear  consistently, 
at  approximate  depths  of  7,  15,  25,  and  41  km.  The  sources  of 
each  of  tlicse  peaks  can  be  identified,  liow'cver,  by  a more  de- 
tailed analysis  of  the  results. 

The  peak  at  7 km,  present  with  an  amplitude  above  the  95s 
significance  level  on  almost  all  of  tlie  depth  plots,  apj'jears 
to  be  a product  of  the  cepstrum  processing.  A peak  corres- 
ponding to  this  delay  time  (approximately  1.8  second)  has 
been  present  on  most  of  the  previous  events,  althougli  usually 
not  witli  such  high  amplitude,  and  is  independent  of  the  cep- 
strum computation  window  length  used.  This  feature  apparently 
arises  from  the  truncation  of  the  central  peak  of  each  cepstrum 
at  a time  of  1.6  second;  techniques  are  being  investigated  for 
removing  the  central  peak  of  the  raw  cepstrums  witliout  intro- 
ducing this  spurious  peak  into  the  depth  plots. 
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Of  tlie  remaining  three  depth  plot  peaks,  the  one  at  25  km 
occurs  most  frequently  above  tlie  95"i  level,  and  provides  the 
best  estimate  of  the  source  depth.  This  conclusion  is  supported 
by  the  observation  of  possible  pi’  arrivals  at  a delay  time  of 
about  7 seconds  on  many  of  the  seismograms. 

The  peak  at  41  km  also  appears  to  represent  a real  seismo- 
gram feature,  since  several  seismograms  show  possible  arrivals 
at  about  10  seconds  after  F.  A 41  km  peak  represents  a delay 
time  of  about  11.3  seconds,  which  is  1 second  longer  than  the 
expected  delay  foi-  sP  from  a 25  km  event.  This  peak  is  most 
St  ra  iglit  fo  rwa  rd  1 y e.xj'laincd  as  a sP  arrival  from  a source 
region  witli  abnormally  low  S velocities.  An  alternative  in- 
terpretation is  that  there  are  two  good  reflectors  above  the 
source,  such  as  water  bottom  and  water  surface,  with  a tAvo- 
way  travel  time  of  about  4 seconds  between  the  two  boundaries. 
Then,  the  25  km  peak  would  be  tlic  result  of  an  echo  off  the 
bottom  reflector,  and  the  41  km  peak  would  come  from  an  echo 
off  the  top  reflector. 

Tl\e  13  km  peak,  representing  a delay  time  of  about  4 sec- 
onds, corresponds  to  the  expected  difference  peak  from  the 
25  km  and  41  km  peaks. 

Summarizing,  the  folloA-.'ing  conclusions  about  the  four 
depth  plot  peaks  have  been  reached: 

• The  7 km  peak  is  a product  of  the  processing. 

• The  13  km  peak  is  a difference  peak  arising 
from  the  25  km  and  41  km  peaks. 

• Two  interpretations  for  the  25  km  and  41  km 
peaks  are  possible: 

1.  The  25  km  peak  represents  a p arrival 
(and  thus  the  event  depth)  , and  the 
41  km  peak  represents  a slow  sP  arrival. 
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2.  Tliere  arc  two  reflectors  above  the 
source,  with  the  2S  km  peak  coming 
from  an  eclio  off  the  deeper  one, 
and  the  41  km  peak  coming  from  an 
echo  off  the  shallower  one. 

This  event  is  a good  one  for  illustrating  liow  many  factors 
can  be  involved  in  obtaining  a depth  estimate.  Seismogram 
plots,  delay  time  information,  and  dejith  plots  from  several 
different  station  groups  were  all  used  in  tlie  interpretation 
of  this  event.  All  this  inrormation  should  be  easily  avail- 
able to  an  analyst  doing  routine  depth  determinations,  and 
each  of  these  features  is  planned  to  be  included  in  future 
program  implementations. 

Remova  1 of  Source  S]iectrum  Structure 

Many  of  the  events  tliat  have  been  processed  through  the 
source  deptli  determination  prograiu  have  given  depth  plots  with 
spurious  peaks  at  shallow  (<15  km)  depths.  These  peaks  could 
be  results  of  structure  in  the  source  wavelet  spectrum,  and 
thus  w'ould  be  eliminated  if  power  spectra  were  inverse  filtered 
with  the  source  wavelet  before  cepstrum  computation.  The 
depth  determination  program  was  modified  to  do  this  inverse 
filtering,  and  two  methods  for  estimating  the  source  spectrum 
were  tried.  The  first  of  these  metliods  averaged  the  spectra 
from  each  cepstrum  computation  window  to  obtain  one  source 
spectrum  per  station,  and  the  second  method  calculated  a sep- 
arate estimate  for  each  window  by  smoothing  the  actual  spectrum 
for  that  window. 

These  technicpies  were  tested  on  both  real  and  synthetic 
data.  On  synthetic  data,  neither  metliod  worked  as  well  as 
conventional  processing,  and,  on  real  data,  neither  method 
worked  at  all.  This  negative  result  is  due  to  the  inadequacy 
of  the  source  spectrum  estimates;  if  a technique  for  obtaining 
accurate  source  spectrum  estimates  could  be  found,  inverse 
filtering  would  probably  have  the  desired  effect. 
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In  order  to  provide  easier  access  to  Seismic  Research 
Observatory  data,  the  source  depth  determination  pro(’,ram  is 
being  implemented  on  the  IBM  TS/44  system  at  tlic  Seismic 
Data  Analysis  tienter.  Since  this  implementation  already  in- 
volves considerable  program  modification,  it  was  decided  to 
rewrite  mud',  of  the  program  as  a ]'>rel imi nary  version  of  the 
analysis  package  that  will  eventually  be  used  for  routine 
source  depth  determination  liy  an  analyst.  Currently,  system 
f 1 ow'cha  r t i ng  has  been  completed  and  coding  lias  been  started 
en  this  new  program  version. 

Tlie  most  significant  difference  between  the  SDAC  imple- 
mentation of  tlie  deptli  determination  ]trogram  and  previous 
versions  is  the  increased  flexibility  for  specification  of 
station  grou]-)ings,  seismogram  time  limits,  and  output  disjilays 
availal-)le  to  the  analyst  in  the  new  SDAC  system.  Tlie  general 
system  flowcliart  (figure  19)  illustrates  how  tliis  is  accom- 
plished. Basically,  the  depth  determination  system  consists 
of  two  steps;  one  that  computes  all  tlie  necessary  cepstrums 
from  tlie  seismogram  input,  and  one  tliat  iiicks  these  cepstrums 
to  generate  depth  jilots.  Cepstrums  cominited  in  tlie  first 
step  are  written  on  the  ceiistrum  file,  which  is  re.ad  by  the 
second  step.  The  second  step  may  also  be  run  as  a separate 
program;  after  viewing  the  first  batch  of  depth  plots  gener- 
ated by  the  system,  the  analyst  may  go  back  to  the  same  cep- 
strum  file  and  generate  new  dejith  plots  using  any  desired 
combination  of  stations,  window  length,  and  total  data  length. 
It  is  planned  that  tliis  generation  of  new  depth  plots  will 
eventually  be  done  on  an  interactive  basis,  without  actually 
submitting  a new  program  run. 

Several  otlier  new  features  have  also  been  incorporated 
into  the  SDAC  depth  determination  system.  These  include  op- 
tional plots  of  individual  cepstrums,  and  optional  cepstrum 
editing  based  on  the  phase  of  each  cepstrum  point. 
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